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1. Introduction 

In this talk, we present the main idea in [|l|] and provide a sampling of the numerical results. 

The Wilson loop operator, IV, is a unitary operator for SU(N) gauge theories and can be used 
as a probe of the transition from strong coupling to weak coupling. Large (area) Wilson loops are 
non-perturbative and correspond to strong coupling. Small (area) Wilson loops are perturbative 
and correspond to weak coupling. 

The probe is defined as 

WN{z,b,L) = {det{z-W)) (1.1) 

and is the characteristic polynomial associated with the operator. W is the Wilson loop operator, 
z is a complex number, N is the number of colors, b = is the lattice gauge coupling and L is 
the linear size of the square loop. (• • •) is the average over all gauge fields with the standard gauge 
action. 

The eigenvalues of W are gauge invariant and so is the characteristic polynomial. The eigen- 
values lie on the unit circle and all of them will be close to unity for small loops. The eigenvalues 
will spread uniformly over the unit circle for large loops. The characteristic polynomial exhibits a 
transition at N = oo when L — > Lc{b). This is a physical transition since Lc{b) will scale properly 
with the coupling, b, as one approaches the continuum limit. 

The scaling function in the double scaling limit can be derived for two dimensional large N 
QCD We have numerically shown that three dimensional QCD falls into the same universality 
class 



2. Two dimensional QCD and a multiplicative matrix model 

Two dimensional gauge theory on an infinite lattice can be gauge fixed so that the only vari- 
ables are the individual plaquettes and these will be independently and identically distributed. 
W = Y\!]=i where Ujs are the transporters around the individual plaquettes that make up the 
loop and « = is equal to the area of the loop. The measure associated with Uj can be set to 

N '-p , tj2 ■ Tj 

P{Uj) = jYe^~ ' ' where IJ j = e ' and £ plays the role of gauge coupling. The dimensionless 
area is given by ? = £^?i which is kept fixed as one takes the limit n — > oo and £ ^ 0. This is called 
the multiplicative matrix model [^. In the continuum limit, the parameters b and L get replaced by 
one parameter, which is denoted by t in the model, and the characteristic polynomial becomes 



3. Average characteristic polynomial 

Using a fermionic representation of the determinant, one can perform the integration over Uj. 
One can then perform the integration over the fermionic variables to obtain the following result for 
the characteristic polynomial: 



QN{z,t) 



SU{N) 
U[N) 



(3.1) 
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Integrating out v gives exact continuum polynomial expressions, 

QN{z,t) = < 



N 1 




k=0 


yk 
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fN 


k=0 


\ k 



SU{N) 



z^-^{-\)^e — U{N) 



(3.2) 



T = t[l + l\ (3.3) 



4. Heat-kernel measure 

The result for QN{z,t) is consistent with the heat-kernel measure for W: 

P{W,T:)dW = Y,dRZR(W)e-'"^'^'^^dW. (4.1) 

R 

R denotes the representation, cIr is the dimension of the representation R and C2{R) is the second 
order Casimir in the representation R. To see this, note that 

QN{z,t) = ifliz-e''^)) = l^z^-'{-\fMt{t). (4.2) 

j=l k=0 

If we now take the average, (• • •) over the heat-kernel measure, we get 

M,(0 = ( X +'^.+-+«.)) = =4.-^^« = (^)e-'^ (4.3) 

i<;i<y2<;3....<A<A' 

5. Zeros of QNiz,t) 

Since W is a unitary operator, the zeros of det(z — W) will lie in the unit circle. One can show 
this remains true for QN{z,t) when the gauge group is SU {N). To see this, we rewrite QN{z,t) for 
SU{N) as 

This is the partition function for a spin model with a ferromagnetic interaction for positive T. 
ln{—z) is a complex external magnetic field. Therefore, the conditions for Lee- Yang theorem [Q] 
are fulfilled and all roots of QN{z,t) lie on the unit circle for SU(N). This is not the case for U(N). 

6. Weak coupling vs strong coupling 

The transition from weak coupling to strong coupling can be intuitively seen using the charac- 
teristic polynomial, 2/v(z,f)- the weak coupling (small area) limit we have t = and QNiz,t) = 
(z — 1)^. Therefore, all roots are at z = 1 on the unit circle. In the strong coupling (large area) limit 
we have t = oo and QN{z,t) =z^ + (—1)^- Therefore, all roots are uniformly distributed on the unit 
circle. 

QN{z,t) is analytic in z for all t at finite N. But, this is not the case as N °° and this leads to 
a transition from weak to strong coupling in the N ^ oo limit. 
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7. Phase transition in an observable - Durhuus-Olesen transition 

There is a critical area, t = 4, such that the distribution of zeros of Qco{z.,t) on the unit circle 
has a gap around z = — 1 for ? < 4 and has no gap for ? > 4 |^ 0. To see this, we note that the 
integral representation (3J) is dominated by the saddle point, v = A(f,z), given by 

ze'('^+2) - 1 

With z = e'^ and w = 2A + 1, p(0) = -g^Re w gives the distribution of the eigenvalues of W on 
the unit circle. 

The saddle point equation at z = — 1 is 

t 

w = tanh-w (7.2) 
showing that w admits non-zero real solutions for t > 4. 



8. Double scaling limit 

AsN ^o°, one can define a scaling region around t = 4 and z = — 1 by 




(8.1) 



a and ^ are the scaling variables that blow up the region near t = 4 and z = — 1 . We can show that 

Mm ' (-ife'-^i-zy'^QMizj) = fjue-"'-""'+^" ^ (8.2) 

which is the scaling function in the double scaling limit associated with the characteristic polyno- 
mial. 

We hypothesize that this behavior in the double scaling limit derived for two dimensional large 
QCD is universal and should be seen in the large N hmit of 3D QCD, 4D QCD, 2D PCM and 
other related models. The modified Airy function, (^{^ , a), is the universal scaling function. 



9. Large N universality hypothesis 

We can now precisely state the continuum large N universality hypothesis that can be numeri- 
cally tested in relevant models. 

Let ^ be a closed non-intersecting loop: {s),s ^ [Oil]- Let 'if (m) be a whole family of loops 
obtained by dilation: x^{s,m) = ^x^(i'),with m > 0. Let W {m,'io{*)) = W {'io{m)) be the family of 
operators associated with the family of loops denoted by '^(*) where m labels one member in the 
family. Define 

Oiv(3;,m,^(*)) = (det(e5 + e'iw {m,'<^{*))) (9.1) 

Then our hypothesis is 



1+ « 



3A^a2('^(*)) 



(9.2) 
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Figure 1: Extrapolation to the continuum limit of the critical area. 



10. Numerical test of the universality hypothesis - 3D large N QCD 

We use standard Wilson gauge action. The lattice coupling b = has dimensions of length. 
We use square Wilson loops of linear length L. We change b to generate a family of square loops 
labeled by L. While doing this, we need to keep 0.42 < b < b{V) where V is the lattice volume 
assumed large enough for large N continuum reduction to hold in the confined phase [||, |6|]. We 
use smeared links in the construction of the Wilson loop operator to avoid corner and perimeter 
divergences. 

We obtain bc{L), ai{L) and a2{L) such that 



lim ^{b,N)ON y 



3N^ J ai{L) 



■bAL) 



1 + 



a 



3Na2{L) 



(10.1) 



This is done by fixing A'^ and L and obtaining estimates for bc{L,N), a\ {L,N) and a2{L,N). We then 
take the limit as N and we check that bc{L), a\ (L) and a2{L) have proper continuum limits as 
L — > oo. The extrapolation to the continuum limit are shown in Fig. [I| Fig. § and Fig. |[ 
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Figure 2: Extrapolation to the continuum limit of the parameter matching to the scaled variable, a. 
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